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Magnesium dihydride (MgH2) stores 7.7 weight % hydrogen, but it suffers from a high thermody- 
namic stability and slow (de)hydrogenation kinetics. Alloying Mg with lightweight transition metals 
(TM = Sc, Ti, V, Or) aims at improving the thermodynamic and kinetic properties. We study the 
structure and stability of Mgi:TMi_2,H2 compounds, x = [0-1], by first-principles calculations at 
the level of density functional theory. We find that the experimentally observed sharp decrease in 
hydrogenation rates for x > 0.8 correlates with a phase transition of Mga;TMi_i,H2 from a fluorite 
to a rutile phase. The stability of these compounds decreases along the series Sc, Ti, V, Cr. Varying 
the transition metal (TM) and the composition x, the formation enthalpy of Mga;TMi_i,H2 can 
be tuned over the substantial range — 2 eV/f.u. Assuming however that the alloy Mg^^TMi-j, 
does not decompose upon dehydrogenation, the enthalpy associated with reversible hydrogenation 
of compounds with a high magnesium content {x = 0.75) is close to that of pure Mg. 

PACS numbers: 71.20.Be, 71.15.Nc, 61.66.Dk, 61.50.Lt 



I. INTRODUCTION 

Hydrogen is a clean energy carrier and an alterna- 
tive to carbon based fuels in the long run;! Mobile 
applications require a compact, dense and safe stor- 
age of hydrogen with a high-rate loading and unload- 
ing capabilityj2i2, Lightweight metal hydrides could sat- 
isfy these requirementsj^i^ Metal hydrides are formed by 
binding hydrogen atoms in the crystal lattice, resulting 
in very high volumetric densities. Reasonable hydrogen 
gravimetric densities in metal hydrides can be achieved 
if lightweight metals are used. 

MgH2 has been studied intensively since it has a rel- 
atively high hydrogen gravimetric density of 7.7 wt. %. 
Bottlenecks in the application of MgH2 are its thermo- 
dynamic stability and slow (de) hydrogenation kinetics. 
These lead to excessively high operating temperatures 
(573 — 673 K) for hydrogen releasei^i^ The hydrogen 
(de)sorption rates can be improved by decreasing the 
particle size down to nanoscalesi ^'^°i^^ It is predicted 
that particles smaller than 1 nm have a markedly de- 
creased hydrogen desorption enthalpy, which would lower 
the operating temperaturejJ^ The production of such 
small particles is nontrivial, however, and the hydrogen 
(de)sorption rates of larger nanoparticles are still too low. 

An additional way of improving the (de)hydrogenation 
kinetics of MgH2 is to add transition metals 
(TMs),^ii^iMii5 Usually only a few wt. % is added, 
since TMs are thought to act as catalysts for the 
dissociation of hydrogen molecules. Recently how- 
ever, Notten and co-workers have shown that the 
(de)hydrogenation kinetics is markedly improved by 
adding more TM and making alloys Mg^-TMji.^;), 
TM=Sc, Ti, X < o.8^^'i^^'20^>22^^^^ The 
basic ansatz is that the rutile crystal structure of MgH2 



enforces an unfavorably slow diffusion of hydrogen 
atoms<21 ScH2 and TiH2 have a fluorite structure, which 
would be more favorable for fast hydrogen kinetics. By 
adding a sufhciently large fraction of these TMs one 
could force the Mg^TM(i_2.) compound to adopt the 
fluorite structure. 

In this paper we examine the structure and stability 
of Mg^TM(i_^)H2, TM = Sc, Ti, V, Cr, compounds by 
first-principles calculations. In particular, we study the 
relative stability of the rutile versus the fluorite struc- 
tures. This paper is organized as follows. In Sec. |TT] 
we discuss the computational details. The calculations 
are benchmarked on the TMH2 simple hydrides. The 
structure and formation enthalpies of the compounds 
Mg^TM(i_a;)H2 are studied in Sec. IIIII and an analysis 
of the electronic structure is given. We discuss the hy- 
drogenation enthalpy of the compounds in Sec. IIVI and 
summarize our main results in Sec. [V] 



II. COMPUTATIONAL METHODS AND TEST 
CALCULATIONS 

We perform first-principles calculations at the level of 
density functional theory (DFT) with the PW91 func- 
tional as the generalized gradient approximation (GGA) 
to exchange and correlation. As transition metals have 
partially filled 3d shells we include spin polarization 
and study ferromagnetic and simple antiferromagnetic 
orderings where appropriate. A plane wave basis set 
and the projector augmented wave (PAW) formalism are 
used,^^'^" as implemented in the VASP code^^i^ The 
cutoff kinetic energy for the plane waves is set at 650 
eV. Standard frozen core potentials are applied for all 
the elements, except for Sc, where we include 3s and 
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FIG. 1: (Color online) (Left) rutile crystal structure of a- 
MgH2, and (right) fluorite crystal structure of q:-TMH2. The 
white spheres represent the hydrogen atoms. 

3p as valence shells, in addition to the usual 4s and 3d 
shells. The Brillouin zone (BZ) is integrated using a reg- 
ular k-point mesh with a spacing ~ 0.02 and the 
Methfessel-Paxton scheme with a smearing parameter of 
0.1 eV.^"^ The self-consistency convergence criterion for 
the energy difference between two consecutive electronic 
steps is set to 10~^ eV. Structural optimization is as- 
sumed to be complete when the total force acting on 
each atom is smaller than 0.01 eV/A. The volumes of 
the unit cells are relaxed, and, where appropriate, also 
their shapes. Finally, we calculate accurate total energies 
for the optimized geometries using the linear tetrahedron 
method,^ 

To calculate the formation enthalpies of the metal hy- 
drides we consider the following reaction. 

xMg+(l-a:)TM + H2(.9) ^Mg,TM(i_,)H2. (1) 

The formation enthalpies (at T = 0) are then obtained 
by subtracting the total energies of the reactants from 
that of the product. A cubic box of size 10 A is applied 
for the H2 molecule. The calculated H— H bond length, 
binding energy, and vibrational frequency are 0.748 A, 
—4.56 eV and 4351 cm~^, respectively, in good agreement 
with the experimental values of 0.741 A, —4.48 eV and 
4401 cm-i,^i^ 

Since hydrogen is a light element, the zero point en- 
ergy (ZPE) due to its quantum motion, is not negligible. 
We find that the correction to the reaction enthalpies of 
Eq. ^ resulting from the ZPEs, is 0.15 ±0.05 eV/Hz, as 
function of the composition x and the transition metal 
TM. Since we are mainly interested in relative formation 
enthalpies, we omit the ZPE energy correction in the fol- 
lowing. 

Before discussing the Mg^TM(i_2,)H2 compounds 
we benchmark our calculations on the simple com- 
pounds MgH2 and TMH2. Under standard conditions 
magnesium- hydride has the rutile structure, a-MgH2, see 
Fig. [2 with space group P42/mnm (136) and Mg and H 
atoms in the 2a and 4/ {x = 0.304) Wyckoff positions, 
respectively. Each Mg atom is coordinated octahedrally 
by H atoms, with two Mg-H distances of 1.94 A and four 
distances of 1.95 A. First row early transition metal hy- 
drides crystallize in the fluorite structure, a-TMH2, see 
Fig. [U with space group Fm3m (225) and TM atoms in 



TABLE L Optimized cell parameters a (c), and calculated for- 
mation enthalpies Ef, of elemental dihydrides in their most 
stable (a) forms. All TMH2 have a fluorite structure, space 
group Fm3m (225), whereas MgH2 has a rutile structure, 
space group P42/mnm (136). 



Compound 


a 

Calc 


[c)A 

Exp 


Ef (eV/f.u.) 
Calc Exp 


MgH2 


4.494 (3.005) 


4.501 (3.010)" 
4.78' 


-0.66 -0.76 


ScH2 


4.775 


-2.09 -2.08 


TiH2 


4.424 


4.454" 


-1.47 -1.45 


VH2 


4.210 


4.27'* 


-0.65 -0.79 


CrHa 


4.140 


3.861'* 





"Ref. 113 
''Ref. [si 
'^Ref. m 
''Rei. EO 

■^Antiferromagneticly ordered. 



a and H atoms in 8c Wyckoff positions. Each TM has 
a cubic surrounding of H atoms with calculated TM-H 
bond lengths of 2.07, 1.92, 1.82 and 1.79 A for Sc, Ti, V, 
and Cr, respectively. By breaking the cubic symmetry by 
hand and reoptimizing the geometry we have confirmed 
that the fluorite structure indeed represents a stable min- 
imum. 

The optimized cell parameters and the calculated for- 
mation enthalpies of the simple hydrides are given in Ta- 
ble m The structural parameters are in good agreement 
both with available experimental data and with previous 
DFT calculations. The formation enthalpies of MgH2 
and VH2 are somewhat underestimated by the calcula- 
tions, whereas those of ScH2 and TiH2 are in excellent 
agreement with experiment. CrH2 is predicted to be un- 
stable with respect to decomposition. 



III. RESULTS Mg^TM(i_^)H2 
A. Structures and formation enthalpies 

Mg^TM(i_^)H2 has the fluorite structure for 2; = 0, 
and the rutile structure for x = 1. We want to establish 
which of the two structures is most stable at intermediate 
compositions x. First we summarize the current status 
of the experimental work on Mgj.TM(i_3.) alloys. 

Experimentally it has been demonstrated that 
Mg^Sc(i_2,) alloys can be reversibly hydrogenated, both 
in thin films, as well as in bulk formj ^^i^''''^^i^°''^^i'^^ 
Mg and Ti do not form a stable bulk alloy, but thin 
films of Mg^Ti(i_3,) have been made, which are read- 
ily and reversibly hydrogenated i^^'^^1^^'^'^1^^'^^1''^^ Thin 
films of Mg^V(i_2;) and Mg^Cr(i_a,) can also be easily 
hydrogenated,— Attempts to produce non-equilibrium 
bulk Mg^Ti(i_2.) alloys by ball milling of Mg and Ti or 
their hydrides have had a limited success so far^ii^i^i^ 
However, MgyTiHy crystals have been made using a 
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high pressure anvil technique42. The same technique 
has been appHed to produce the hydrides MggVHj^ and 
MgaCrHj^.i^i^ia. 

The crystal structure of Mg^Sc(i_2.)Hy and 
Mg^Ti(i_j.)Hj, in thin films, x < 0.8, y ~ 1-2, is 
cubic, with the Mg and TM atoms at fee positions. No 
detectable regular ordering of Mg and TM atoms at 
these positions has been foundj^i^i^ In contrast, the 
Mg and TM atoms form simple ordered structures in the 
high pressure phases.— '^'^'^'^ The hydrogen atoms in 
Mgo.65Sco.35Hj,, y « 1-2, assume tetrahedral interstitial 
positions, as is expected for the fluorite structure^ii^ 
In the MgyTiHie high pressure phase the metal atoms 
are in fee positions and are ordered as in the Ca7Ge 
structure42. The H atoms are in interstitial sites, but 
displaced from their ideal tetrahedral positions. 

The latter structure can be used as a starting 
point to construct simple, fluorite-type structures for 
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p-MgH, 



Mg,TM(i_,)H2, < X < 1. 
use the C&rGe structure to order the metal atoms, for 
X = 0.25, 0.75 the CuaAu (LI2) structure, and for x = 0.5 
the CuAu (ilo) structure. The H atoms are placed at 
or close to tetrahedral interstitial positions. To model 
Mg^TM(-]^_^)H2 in rutile-type structures we use the a- 
MgH2 structure as a starting point. For each compo- 
sition X we choose the smallest supercell of that struc- 
ture corresponding to that composition. As the atomic 
volumes of the various TM atoms differ, in each of the 
structures and compositions the cell parameters are op- 
timized, as well as the positions of all atoms within the 
cell. Care is taken to allow for breaking the symmetry in 
the atomic positions. In particular the hydrogen atoms 
are often displaced from their ideal tetrahedral positions. 
Although these structures are then no longer ideal fluo- 
rite structures anymore, we still use the term fluorite in 
the following. Test calculations show that the total ener- 
gies of the hydrides of (quasi-)random alloys are similar 
to those of the simple ordered structures. In particular, 
the relative stability of the fluorite vs. the rutile struc- 
tures is not extremely sensitive to the relative ordering 
of the metal atoms. 

The calculated volumes T^(Mg^TM(i_3;)H2) of the flu- 
orite structures, normalized per formula unit (f.u.), are 
shown in Fig. [5) The volumes of a few bulk compounds 
and compositions can be extracted from experimental 
data, thus providing a check on the calculations. Inter- 
polating the results in Ref . l4ll gives ^(MgQ 55800.35112) — 
27.5 A'^, whereas the in terp olated calculated value from 
Fig. His 26.9 A3. Ref. Ha g 

ives l^(Mgg g7gTio.i25 Ha) ^ 
27.3 A'^, compared to the calculated value 26.6 A''. 
These differences between experimental and calculated 
volumes are consistent with the differences between ex- 
perimental and calculated lattice parameters of the sim- 
ple hydrides, see Table HI Experiments on thin fllms 
give y(MgQ 7Tio.3H2) — 26.4 A^^, if one assumes cu- 
bic symmetryj21 The interpolated computational value 
is 24.4 A"^. The difference between these values is not 
excessively large, but considering the smaller difference 



Pqj. 2; = 125 875 we FIG. 2: (Color online) The volumes per formula unit in of 



Mg^TM(i_j;)H2 in the fluorite structure, as a function of the 
composition x for TM = Sc, Ti, V, Cr (from top to bottom). 



found for the bulk composition Mgp 875Tio.i25H2, it may 
suggest that the structure of thin films is slightly different 
from that of bulk. 

The trends observed in Fig. [2] can be interpreted 
straightforwardly. The cell volumes of a-ScH2 and the 
cubic /3-MgH2 structure^ are within 1.4% of one an- 
other, which explains why the volumes calculated for 
Mgj.Sc(i_j;)H2 only weakly depend on the composition 
X. The cell volumes of the other TMH2 are smaller, 
hence one expects the volumes I/(Mg^TM(i_^)H2) to 
increase with x. At fixed composition x, the volumes 
y(Mgj.TM(i_^)H2) decrease along the series Sc, Ti, V, 
Cr, as the atomic volumes of the TMs decrease corre- 
spondingly. According to Zen's law of additive volumes 
one would expect^ 

l/(Mg,TM(i_,)H2) = x!/(/3-MgH2) + (1 - x)V{TMR2). 

(2) 

The curves shown in Fig.[2]deviate slightly, but distinctly, 
from straight lines, with a maximum deviation of 5 %. 
This deviation is consistent with the experimental ob- 



servations on Mg^Ti(i_2.)H2i^'^ It is also observed in 
simple metal alloysi^ 

The calculated formation enthalpies of Mg2,TM(i_j.)H2 
are shown in Fig. [3l Clearly for all TMs the fluorite 
structure is more stable than the rutile structure for all 
X smaller than a critical value, Xc- The critical composi- 
tion Xc at which the rutile structure becomes more sta- 
ble can be guessed by interpolation and is in the range 
Xc « 0.8-0.85 for all TMs. The fact that the critical 
composition is fairly high might be guessed from the en- 
ergies of the fluorite and rutile structures of the simple 
hydrides. As a first estimate of the critical composition 
Xc below which the fluorite structure is stable one may 
try a linear interpolation between the pure TMH2 com- 
pounds, X = 0, and MgH2, x = 1. The rutile structure 
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FIG. 3: (Color online) The formation enthalpy (per formula 
unit) of the Mg^TM(i_^)H2 compounds as obtained from spin 
polarized calculations. The values for the fluorite and rutile 
structures are represented by squares (solid lines) and trian- 
gles (dashed lines), respectively. 



of TMH2 is more unstable than the fluorite structure 
by A(TMH2) = 0.75, 0.65, 0.68, 0.36 eV/f.u. for TM 
— Sc, Ti, V, Cr, respectively. The difference in forma- 
tion enthalpy between the a (rutile) and /? (cubic) phases 
of MgH2 is A(MgH2) = 0.10 eV/f.u. Linear interpola- 
tion then gives = A(TMH2)/(A(TMH2) A(MgH2)), 
which results in Xc ~ 0.9 for Sc, Ti and V, and x^ — 0.8 
for Cr. Whereas these values may seem a good first guess, 
they are somewhat too high as compared to the crossing 
points Xc observed in Fig. [3l Moreover, as this figure 
shows, in particular the curves for Ti and V are far from 
linear, so the results for the linear interpolation may be 
somewhat fortuitous. 

Experimental results indicate that Mgo.7Tio.3H2 has 
the fluorite structure and Mgo.9Tio.1H2 has the rutile 
structurei^^ This agrees with the results shown in Fig. [U 
where a fluorite to rutile phase transition takes place at 
the composition x = Xc 0.83. It has been suggested 
that the fluorite structure allows for a much faster kinet- 
ics of hydrogen loading and unloading. Experimentally 
it has been observed that the dehydrogenation kinetics of 
Mg^Sc(i_j;) and Mg^Ti(i_3.) becomes markedly slower if 
X > xq — 0.8 . The results shown in Fig. [T] suggest that 
xo = Xc, i.e. the composition at which the phase transi- 
tion between fluorite and rutile structures takes place. 

The dehydrogenation kinetics of Mg2,Ti(i_2,) also be- 
comes gradually slower with decreasing x, for x < Xc, i.e. 
where the compound remains in the fluorite structure, 
although it is still faster than for x > Xc^^i^i Whereas 
kinetic studies are beyond the scope of the present pa- 
per, we speculate that a volume effect might play a role 
here. In the fluorite structure the hydrogen atoms oc- 
cupy interstitial positions close to the tetrahedral sites. 



Diffusion of hydrogen atoms is likely to take place via 
other interstitial sites such as the octahedral sites. The 
smaller the volume, the shorter the distance between 
such sites and the occupied positions, or in other words, 
the shorter the distance between a diffusing hydrogen 
atom and other hydrogen atoms in the lattice. This 
may increase the barrier for diffusion. As the volume 
of Mg^V(i_a;) and Mg^Cr(i_2,) is smaller than that of 
Mgj.Sc(i_j;) and Mg^Ti(i_2,) (at the same composition 
x), this might also explain why the dehydrogenation ki- 
netics of the former compounds is much slower Jii We note 
that the smaller volume of Mg^V(i_a;) and Mgj.Cr(i_a;) 
is accompanied by a distortion of the structures consis- 
tent with the limited space available to accommodate 
the hydrogen atoms. For instance, in Mgp 75V0.25H2 and 
^60.75^^^0. 25H2 the hydrogen atoms are displaced con- 
siderably from the tetrahedral positions, and the coordi- 
nation number of V and Cr (by hydrogen) is 7, instead 
of 8 as in case of a perfect fluorite structure. 



B. Electronic structure 

To analyze the electronic structure of the compounds 
Mgj.TM(i„2,)H2, we start with the density of states 
(DOS) of the pure hydrides a-MgH2 and TMH2 as shown 
in Fig. m The bonding in MgH2 is dominantly ionic; oc- 
cupied hydrogen orbitals give the main contribution to 
the valence states, whereas the conduction bands have a 
signiflcant contribution from the Mg orbitalsi^ As usual, 
ionic bonding between main group elements results in an 
insulator with a large band gap. In contrast, the tran- 
sition metal dihydrides are metallic, as demonstrated by 
Fig. m The peak in the DOS at low energy, i.e. be- 
tween —9 and —2 eV in ScH2 to between —12 and — 4 eV 
in CrH2, is dominated by hydrogen states. The broad 
peak around the Fermi level consists of transition metal 
d-states. It suggests that bonding in TMH2 is at least 
partially ionic. The TM s-electrons are transferred to 
the H atoms, whereas the d-electrons largely remain on 
the TM atoms. The DOSs of TMH2, TM=Sc, Ti, V, Cr, 
are very similar in shape. As the number of d-electrons 
increases from one in Sc to four in Cr, the Fermi level 
moves up the d-band in this series. As the DOS at the 
Fermi level increases, it enhances the probability of a 
magnetic instability. Indeed we flnd CrH2 to be antifer- 
romagnetic with a magnetic moment of 1.5/iB on the Cr 
atoms. The antiferromagnetic ordering is 51 meV/f.u. 
more stable than the ferromagnetic ordering, which is 7 
meV/f.u. more stable than the non-polarized solution. 
In the other TMH2 we do not find magnetic effects. 

Additional information on the type of bonding can be 
obtained from a Bader charge analysis. In Q;-MgH2 the 
Bader charges are Qug = +1.59e (and Qh = — 0.80e, 
since the compound is neutral), which confirms that this 
compound is dominantly ionic. The results for TMH2 
are shown in Table HTl They indicate that the ionicity in 
ScH2 is comparable to that in MgH2. Furthermore, the 
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FIG. 4: (Color online) Densities of states of MgH2 and TMH2 
(left column) and of Mgy 75TM0.25H2 (right column) for TM 
= Sc, Ti , V, Cr. The shaded areas give the projected densities 
of states on the TM d orbitals. For CrH2 the nonmagnetic 
DOS is given for simplicity reasons; CrH2 is antiferromagnetic 
(see text). 



ionicity decreases along the series Sc, Ti, V and CrH2. 
Comparison to Table[T]shows that the decrease in ionicity 
correlates with a decrease in formation enthalpy. 

These results on the simple hydrides help us to analyze 
the electronic structure and bonding in Mg^TM(i_2;)H2. 
We show results for the fluorite structure only, since that 
is the more stable structure over most of the composi- 
tion range. As an example. Fig. 2] shows the calculated 
DOSs of Mgo 75TM0.25H2. One can qualitatively inter- 
pret these DOSs as a superposition of the DOSs of MgH2 
and TMH2. The bonding states at low energy, compris- 
ing the first broad peak in the DOS, consist mainly of 
filled hydrogen states. The peaks close to the Fermi level 
are dominated by TM d-states. The Fermi level moves 
up the d-band through the series Sc, Ti, V, Cr. At higher 
energy we find the (unoccupied) Mg s states. The basic 
structure of the DOSs remains the same for all composi- 
tions Mg2,TM(i_a;)H2. As X increases, the TM d contri- 
bution of course decreases. In addition, the TM d-peak 
becomes narrower with increasing x, as the distance be- 
tween the TM atoms increases. 



TABLE 


II: Bader 


charge 


analysis of 


TMH2 and 


00.75 

e. 


0.25H2. All 


charges 


Q are given in units of 




TMH 


2 


Mgo. 


,5TMo.25H2 


TM 


Qtm (e) 


Qh (e) 


Qtm (e 


) Qh (e) 


Sc 


+1.51 


-0.75 


+ 1.57 


-0.80 


Ti 


-fl.17 


-0.59 


+ 1.18 


-0.76 


V 


+1.09 


-0.55 


+0.98 


-0.73 


Cr 


+0.89 


-0.45 


+0.68 


-0.69 



Narrowing of the d peak can give rise to magnetic in- 
stabilities. The tendency to such instabilities increases 
along the series Sc, Ti, V and Cr. The development of 
nonzero magnetic moments of course strongly depends 
upon the structure. Nevertheless, for Mg^Sc(i„j;)H2 and 
Mgj.Ti(]^_j.)H2 we see a tendency to form magnetic mo- 
ments on the TMs only if a; > 0.8. For Mg^V(i_a;)H2 
this occurs if a; > 0.5, and for Mg^Cr(i_2,)H2 one can 
find magnetic instabilities over the whole composition 
range. Most of the structures have a finite DOS at the 
Fermi level, which, might indicate a metallic behavior. 
One cannot conclude this on the basis of a DOS alone, 
however, but should also critically evaluate possible lo- 
calization and on-site correlation effects. There are a few 
exceptions. In particular cases low spin states can be 
more stable, such as for MgQ 75Cro.25H2 in the fiuorite 
structure. Cubic crystal field splitting by the hydrogens 
surrounding the Cr atom results in a gap between Cg and 
t2g states, the eg states being lowest in energy. The lat- 
ter are filled by the four d electrons of Cr, which makes 
this particular structure insulating, see Fig. 01 The DOS 
of Mgg 75V0.25H2 in the fluorite structure is explained 
by the same mechanism. However, as V only has three 
d electrons, each V atom obtains a magnetic moment 
of 1 fj,B- The distance between the TM atoms is fairly 
large in most compositions that have nonzero magnetic 
moments, which suggests a small magnetic coupling be- 
tween the TM atoms, a low Neel or Curie temperature, 
and paramagnetic behavior at room temperature. Excep- 
tions are the Cr compounds with a substantial amount 
of Cr, as discussed above. 

A Bader charge analysis of Mg^TM(i_3.)H2 can be 
made, similar to the simple hydrides. For all compo- 
sitions QMg ~ +1.6e, i.e. close to the value found in 
Q!-MgH2. As an example the Bader charges on the TM 
and H atoms in MgQ 75TM0.25H2 are given in Table [TTl 
The charge on the TM atoms decreases along the series 
Sc, Ti, V, and Cr as in the simple hydrides, but com- 
pared to the latter, it is somewhat smaller on V and 
Cr. The charges on the H atoms in Mgg 75TM0.25H2 are 
roughly the proportional average of the charges on the H 
atoms in MgH2 and TMH2 . The charge analysis of the 
Mgj.TM(i_a;)H2 compounds is consistent with the bond- 
ing picture extracted from the DOSs. 
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IV. DISCUSSION 

We discuss to what extend the Mg-TM aUoys are 
suitable as hydrogen storage materials. The forma- 
tion enthalpies of Mg^TM(i_2,)H2 are shown in Fig. [31 
Lightweight materials require a high content of mag- 
nesium, but to have a stable fluorite structure it 
should not exceed the critical composition Xc, as dis- 
cussed in Sec. IIII Al We focus upon the composition 
MgQ 75TM0.25H2 in the following discussion. The calcu- 
lated formation enthalpies arc —0.83, —0.59, —0.43 and 
-0.30 eV/f.u. for TM = Sc, Ti, V, and Cr, respectively. 
For applications the binding enthalpy of hydrogen in the 
lattice should be < 0.4 eV/H2,^''^''* which indicates that 
the Sc and Ti compounds are too stable. The formation 
enthalpies of the V and Cr compounds could be in the 
right range. However, the parameter that is most rele- 
vant for hydrogen storage is the hydrogenation enthalpy. 
Assuming that the alloy does not dissociate upon dehy- 
drogenation, the hydrogenation enthalpy corresponds the 
reaction 



Mg,TM(i 



H2(g) -^Mg,TM(i_,)H2. (3) 



To assess the hydrogenation enthalpy, one can break 
down the formation enthalpy associated with Eq. ([T]) into 
components, similar to the decomposition used in Ref.ls?!. 
We write the formation enthalpy as a sum of three terms. 

(i) The enthalpy required to make the Mg-TM alloy in the 
fee structure from the elements in their most stable form. 

(ii) The energy required to expand the fee lattice in order 
to incorporate the hydrogen atoms, (iii) The energy as- 
sociated with inserting the hydrogen atoms. The results 
of this decomposition for Mgg 75TM0.25H2 are given in 
Fig. [5l To facilitate the discussion, a similar decomposi- 
tion is shown for the simple hydrides, where (i) only con- 
sists of transforming the pure metal into the fee structure. 
In contrast to Ref. [53, we use the spin-polarized fee alloy 
for calculating the contributions (i) and (ii), as this will 
make the extraction of the hydrogenation enthalpy easier. 
In the cases where the magnetic moment is nonzero, we 
study both ferromagnetic and antiferromagnetic order- 
ing. As for the simple hydrides, Cr compounds generally 
have an antiferromagnetic ordering. 

The lattice expansion energy (ii) of the compounds 
Mgg 75TM0.25H2 is < 0.1 eV for all TMs, see Fig.EJb). 
It is in fact comparable to that of pure Mg, see Fig. [5fa). 
At the composition Mgg 75TM0.25H2, the effect on the 
energy of changing the unit cell volume is dominated 
by Mg. For these compounds the lattice expansion only 
plays a minor role in the formation energy, in contrast to 
the simple hydrides, where the lattice expansion gives a 
significant contribution. The hydrogen insertion energies 

(iii) are also remarkably similar for the Sc, Ti, and V 
compounds. Again this is in sharp contrast to the corre- 
sponding energies for the simple hydrides, which strongly 
depend on the TM. The hydrogen insertion energies for 
the compounds are in fact similar to that of pure Mg. At 
the composition Mgg 75TMg.25H2 also this energy is then 
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FIG. 5: Decomposition of the formation energy into 
(i) tlie formation energy of the spin-polarized fee metal 
'^8o.75TMo.25 (black), (ii) the lattice expansion energy (gray), 
and (iii) the hydrogen insertion energy (white) ; (a) the simple 
hydrides MgHa, TMH2; (b) Mg^ 7gTMo.25H2. 



dominated by Mg. Only the compound Mgg 75Cro.25H2 
has a somewhat smaller hydrogen insertion energy. The 
reason for this is that the energy gained by magnetic 
ordering of the alloy Mgg 75Crg.25 is relatively high, as 
compared to the other compounds. This contribution 
stabilizes the alloy with respect to the hydride, which is 
nonmagnetic. 

The formation enthalpy of the fee alloys Mgg 75TM0.25 
(i) shows the largest variation as a function of the TM, 
relative to the contributions (ii) and (iii). Whereas the 
alloy formation energy is negative for TM = Sc, indicat- 
ing that this alloy is stable, it is positive for Ti, V, and Cr, 
meaning that these alloys are unstable. This result agrees 
with the experimental finding that of the Mg-TM alloys 
considered here, only a stable Mg-Sc alloy exists in bulk 
form. The substantial increase of the alloy formation en- 
thalpy in the series Sc, Ti, V, is largely responsible for 
the variation of the formation energy of the correspond- 
ing hydrides Mgg 75TMg.25H2. The alloy formation en- 
ergy of Mgg 75Crg.25 is similar to that of Mgg 75 Vq. 25, due 
to a relatively high spin-polarization energy, as discussed 
in the previous paragraph. 

The hydrogenation enthalpy according to Eq. ^ can 
be determined by summing the contributions (ii) and (iii) 
of Fig. [5l Since the most stable structure of the alloys is 
not always the fee structure, one should however subtract 
the energy required to convert the alloys from their most 
stable structure to an fee structure. We find, for instance. 
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that for Mgg 75X19.25 the foe structure is 0.04 eV/f.u. less 
stable than the hep structure. Indeed thin film experi- 
ments on Mg^Ti(i_j.) yield yield an hep structure.— 
For Mgg 75SC0.25 the fee structure is more stable than 
the hep structure by 0.05 eV/f.u.. 

The calculated hydrogenation enthalpy of Mgp 75800.25 
is —0.79 eV/f.u., in good agreement with the experimen- 
tal value of —0.81 eV/f.u..^^ The calculated hydrogena- 
tion enthalpy of Mgg 75Tio.25 is —0.76 eV/f.u., which is 
in good agreement with the experimental value of —0.81 
eV/f.u. of Ref. [2^, obtained if the thin film correc- 
tion suggested there is included. These hydrogenation 
enthalpies are remarkably similar to that of pure Mg, 
strongly suggesting that alloying Mg with these TMs does 
not improve this energy as compared to pure Mg. The 
most stable structures of Mgg 75TM0.25, TM = V, Cr, are 
not known, but judging from the Sc and Ti compounds 
the energy difference between the fee and the most sta- 
ble structures will be small. Neglecting this energy dif- 
ference upper bounds for the hydrogenation enthalpies 
of Mgg 75V0.25 and Mgg 75^0.25 are -0.72 and —0.57 
eV/f.u., respectively. Again this indicates that alloying 
Mg with these TMs does not improve the hydrogenation 
enthalpy substantially. 

V. SUMMARY 

In summary, we have studied the structure and sta- 
bility of Mg^TM(i_^)H2, TM = Sc, Ti, V, Cr, com- 
pounds by first- principles calculations. We find that for 



X < Xc ^ 0.8 the fluorite structure is more stable than 
the rutile structure, whereas for x > Xc the rutile struc- 
ture is more stable. This phase transition correlates with 
the observed slowing down of the (de)hydrogenation ki- 
netics in these compounds if x exceeds the critical com- 
position Xc- The density of states of these compounds 
is characterized by the valence bands being dominated 
by contributions from the hydrogen atoms, wherease the 
TMs have partially occupied d states around the Fermi 
level. As x increases and/or one moves down the TM se- 
ries, the tendency to magnetic instabilities increases. The 
formation enthalpy of Mgj;TMi_a;H2 can be tuned over a 
substantial range, i.e. 0-2 eV/f.u., by varying TM and x. 
To a large part this reflects the variation of the formation 
enthalpy of the alloy Mg^TM(i_3,), however. Assuming 
that the alloys do not decompose upon dehydrogenation, 
the hydrogenation enthalpy then shows much less vari- 
ation. For compounds with a high magnesium content 
{x = 0.75) it is close to that of pure Mg. 
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